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Abstract 

The cross section of the process e + e~ — ► 7r°7r°7 has been measured in the cm. energy 
range 600-970 MeV with the CMD-2 detector. The following branching ratios have 
been determined: B{p° -» 7r°7r°7) = (5.2±J;| ± 0.6) x 10" 5 and B(w -» vr vr 7) = 
(6.4±jj;3 ± 0.8) x 10~ 5 . Evidence for the p° -» f (600)7 decay has been obtained: 
B{p° -» /o(600)7) = (6.0j|;? ± 0.9) x 10~ 5 . From a search for the process e+e" -> 
r/7r°7 the following upper limit has been obtained: £>(u — ► r]TT ^/) < 3.3 x 10~ 5 at 
CL. 
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1 Introduction 



Radiative transitions of vector mesons into two pseudoscalar mesons have 
been attracting attention since long ago as a possible test of various low en- 
ergy theoretical models and a source of information on controversial scalar 
states [1,2,3,4,5,6,7,8,9,10,11,12,13]. After the reliable observation of the / (980) 
and a (980) states in the 0(1020) meson decays by SND [14] and CMD-2 [15], 
recently confirmed by KLOE [16], the interest moved to the p and u> meson 
decays. 

Information on the p(u) decays to the 7in'j(r]7i 'y) final states is rather scarce. 
Because of the large background from the initial state radiation in the process 
e + e~ — > 7r + 7r~, a search for p{uj) decays into the 7r + 7r~7 final state is difficult 
and among many such experiments [17] only one succeeded in the observation 
of the decay p° — > 7r + 7r~7 [18]. A long search for the uo — > 7r°7r°7 decay (see [17] 
and references therein) finally proved successful for the GAMS Collaboration, 
which observed it in ir~p collisions with the branching fraction of (7.4 ±2.5) x 
10~ 5 [19]. Recently a high statistics study of the p(oo) energy range has been 
performed by the SND Collaboration [20,21]. They observe both p and u 
decays into 7r°7r°7 with a branching ratio higher than that predicted by vector 
dominance. While for the p meson the excess can be explained by the p° — > 
/o(600)7 decay, first evidence for which is reported by SND, the reasons for 
the higher probability of the corresponding uj decay are not yet clear. 

In our recent paper [22] we described a study of the process e + e _ — > 7r°7r°7 in 
the cm. energy range 920-1380 MeV, i.e. above the threshold of um° produc- 
tion, using the CMD-2 detector at the VEPP-2M e + e~ collider. In this work we 
report on the measurement of the cross section of the process e + e~ — > 7r°7r°7 
in the cm. energy range 600-970 MeV with CMD-2. Also described is the first 
search for the process e + e~ — > rjn®^ in this energy range. 



2 Experiment 

The general purpose detector CMD-2 has been described in detail elsewhere [23] 
Its tracking system consists of a cylindrical drift chamber (DC) and double- 
layer multiwire proportional Z-chamber, both also used for the trigger, and 
both inside a thin (0.38 X ) superconducting solenoid with a field of 1 T. The 
barrel Csl calorimeter with a thickness of 8.1 X placed outside the solenoid 
has energy resolution for photons of about 9% in the energy range from 100 to 
700 MeV. The angular resolution is of the order of 0.02 radians. The end-cap 
BGO calorimeter with a thickness of 13.4 X placed inside the solenoid has 
energy and angular resolution varying from 9% to 4% and from 0.03 to 0.02 
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radians, respectively, for the photon energy in the range 100 to 700 MeV. The 
barrel and end-cap calorimeter systems cover a solid angle of 0.92 x 4n radians. 

This analysis is based on a data sample corresponding to integrated luminosity 
of 7.7 pb" 1 collected in 1998-2000 in the energy range 600-970 MeV. The step 
of the cm. energy scan varied from 0.5 MeV near the u peak to 5 MeV far 
from the resonance. The beam energy spread is about 4 x 10 -4 of the total 
energy. The luminosity is measured using events of Bhabha scattering at large 
angles [24]. 

We use Monte Carlo simulation (MC) to model the response of the detector 
and determine the efficiency. Due to the beam background additional ( "fake" ) 
clasters can appear in the calorimeter. The corresponding probability as well 
as photon energy and angular spectra are obtained directly from the data 
using the process e + e~ — > 7r + 7r~7r°. Then these photons are mixed with the 
detector response during simulation. 



3 Data analysis 



At the initial stage, we select events which have no tracks in the DC, five pho- 
tons, the total energy deposition 1.7 < £tot/-Ebeam < 2.2, the total momentum 
Ptot/E hcam < 0.3 and at least three photons detected in the Csl calorimeter. 
The minimum photon energy is 30 MeV for the Csl and 40 MeV for the BGO 
calorimeter. Then a kinematic fit requiring energy-momentum conservation 
is performed with an additional reconstruction of two n° mesons. We require 
good reconstruction quality (x 2 < 5) and the ratio of the reconstructed to 
measured energy to be 0.75 < uoi / E { < 1.8 for each photon. After this stage 
219 events remain in the whole energy range. 

The dominant background comes from the processes 

e + e~ — > 777, r] — > 37T°, (1) 

e + e- -> 7r°7, (2) 

e + e- -> 37, 4 7 . (3) 



Events from the process (1) can imitate signal events if two soft photons are 
lost. The processes (2,3) can meet the selection criteria with one ore two 
additional ( "fake" ) photons coming from shower splitting, "noisy" electronic 
channels in the calorimeter and beam background. 

To determine the background contribution, the following procedure is used. 
All above listed processes are simulated using the Monte Carlo. To decrease 
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Fig. 1. The E tot (left) and \ 2 distributions. The points with error bars represent ex- 
perimental events, hatched histograms show the MC simulation for the background 
processes and open histograms are a sum of the signal and background MC. The 
arrows indicate the cuts imposed. 



a statistical uncertainty, ten times more events than is expected from the 
background cross section and luminosity are used at each energy. Then the 
same selection criteria as for the data are applied. The obtained number of 
selected events is divided by ten and subtracted from experimental data at 
each energy. In total, about 29 background events are expected after this 
procedure. 

Figure 1 shows the E tot and x 2 distribution for the data, signal and background 
MC. We use the x 2 distributions to estimate the accuracy of the background 
estimation. The experimental distribution was fitted by a sum of MC and 
background. The ratio N^/N^p = 1.2 ±0.2 was obtained. We conclude that 
the background level is estimated well and its systematic error does not exceed 
40%. This results in a 6% systematic uncertainty in the cross section. 



3. 1 Approximation of the cross sections 



At each energy point the cross section of the process a is calculated from 
the observed number of events and background MC expectation using the 
following formula: 

^ = W' (4) 
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E, MeV E, MeV 

Fig. 2. The cross section of the process e + e~ — ► 7r°7r°7. The points with error bars 
represent the experimental data, while the solid and dashed curves correspond to 
the results of the fit I and III, respectively. 

where iV ex p is the number of observed events, A^ g is the expected number of 
background events from MC, L is the integrated luminosity, e is the detection 
efficiency and (1 + 5) is the radiative correction at the corresponding energy. 

To calculate the detection efficiency, we use Monte Carlo simulation taking into 
account the neutral trigger (NT) efficiency. NT is based on the information 
from the Csl calorimeter and its efficiency depends on the number of clusters 
and total energy deposition. The NT efficiency is estimated using events of the 
process e + e~ — > 7r + 7r~7r°. We require the charged trigger signal and three or 
more clusters in the Csl calorimeter, and study the NT efficiency as a function 
of the energy deposition in Csl. The NT efficiency varies from 85% at cm. 
energy 600 MeV to about 95% at 980 MeV. 

The radiative corrections are calculated according to [27]. The dependence of 
the detection efficiency on the energy of the emitted photon is determined 
from simulation. 

The obtained Born cross section of the process e + e~ — > 7r°7r°7 is shown in 
Fig. 2 while Table 1 lists detailed information on the analysis of this reaction. 
No events were selected in the energy range 600 to 690 MeV, therefore our 
results are presented as upper limits at 90% C.L. The Feldman-Cousins proce- 
dure [25] was used to calculate errors (upper limits) at each energy. This cross 
section (the "dressed" one from the column VII) was used in the approxima- 
tion of the energy dependence with resonances. 

Meanwhile, for applications to various dispersion integrals like that for the 
leading order hadronic contribution to the muon anomalous magnetic mo- 
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Table 1 

The cm. energy, integrated luminosity, detection efficiency, number of observed 
events, background expectation, radiative correction, Born cross section cr, vacuum 
polarization correction and "bare" cross section a of the process e + e~ — ► 7r°7r°7. 
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merit, the "bare" cross section should be used. Following the procedure in 
Ref. [26], the latter is obtained from the "dressed" one by multiplying it by 
the vacuum polarization correction |1 — n(s)| 2 , where IT(s) is the photon po- 
larization operator calculated taking into account the effects of both leptonic 
and hadronic vacuum polarization. The values of the correction and the "bare" 
cross section a are presented in two last columns of Table 1. 

The maximum likelihood method is applied to fit the experimental data ob- 
tained from the relation (4). We parameterize the amplitude of the process by 
a sum of the p and u contributions. The former contains the p — > tun tran- 
sition plus one more mechanism beyond the vector dominance model which 
is chosen to be the p° — > fo (600)7 one - Because of the small width of the u 
meson, the p — > un amplitude is rapidly falling below the utt° threshold, 
making thereby the p° — > luit and p° — > f (600)7 transitions distinguishable. 
On the contrary, the small width of the oj meson prevents from distinguishing 
various mechanisms possibly existing in the u — > 7r°7r°7 decay. For this reason 
we parameterize the u meson amplitude by the u — > pn transition only. 

The Born cross section of the process is written as: 

0-7T°7r° 7 ( S ) = / |ArO«o 7 ( s )| 2 ^ » ( 5 ) 



where d$> is the final state phase space and 

2 2 2 9 

777, m m m 

A 7r 7r 7 = A p ^ UJ7r o(—^- + + A>-»/ (600)7 7p + ^^7r ( '7r" 7 TP ■ (6) 

Up U pi Up U w 



Here the first term describes the amplitude of the e + e~ — > p, pi — > ujtx tran- 
sition, while the second and third ones are the e + e~ — > p — > / (600)7 and 
e + e~ — > — > p7T° amplitudes, my is the mass and Dy is the propagator of 
the vector meson V given by D v (s) = s — m v + iy/sT v (s), T v (s) is the corre- 
sponding energy dependent width. The real parameter a = g p!u)ll l 1 g pijm is the 
ratio of the coupling constants for the p and pi mesons. The A p _ tU)7r o ampli- 
tude, proportional to the coupling constant g pu}n of the p — > uj-k transition, is 
written as in our previous analysis above 1 GeV [22]. 

The coupling constant g pum and the following branching ratios are used as fit 
parameters: 

B(p° - /o(600) 7 -> 7rV 7 ) = - / |^ /o(600h (m p )| 2 rf$ , (7) 

Up J 

B{u -> 7T 7r°7) = — f |A UJ ^ 7r o 7r o 7 (m w )| 2 rf$ . 
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Table 2 

The fit results in various models 
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Then the total branching fraction of the p° — > 7r°7r°7 decay is calculated from 
the following formula: 



B(p -> 7T 7r°7) = — j \A p ^o(m p ) + A p ^ Mmoh (m p )\ 2 d$ . (8) 



In (7) and (8) cry is the cross section at the resonance peak without taking 
into account other contributions: 

, 2 , 12ttB{V -> e + e~)B{V -> 7r°7r°7) 

CTy = CTe+e-^y^o^^myJ = j ' ( 9 ) 



— > e + e ) and — > 7r°7r°7) are the corresponding branching ratios. 



3.2 Results of the fits 



In all the following fits the p, 00 and pi meson masses and widths are fixed at the 
world average values [17]. The / (600) mass and width are badly known [17]. 
Therefore, we use a wide range of these parameters: Mj ( 600 ) = 400-800 MeV, 
r/ (6oo) = 300-600 MeV. For the p and pi resonances the energy dependence 
of the total width was described similarly to Ref. [22] while for the u meson 
the total width was assumed to be energy independent. We perform three 
main fits: with g puJ7T equal to the value (16.7 ± 0.4 ± 0.6) GeV -1 obtained in 
our analysis above 1 GeV [22] (fit I), with free g pu)7T (fit II) and without a 
contribution from the p° — * f (600)7 decay (fit III). The results of the fits are 
shown in Table 2 and in Fig. 2 by the curves. 

The value of the coupling constant g puJ7r = 16.2 ±1.4 obtained in the second fit 
is in good agreement with the one from our measurement of the e + e~ — > ujtt 
cross section above 1 GeV [22]. 

The fits taking into account the p° — > f (600)7 decay mode (fits I and II) better 
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describe data. The branching fraction B(p° — > /o (600)7) differs from zero 
by two standard deviations. However, the fit III also has good x 2 / 11 - d. f- = 
26.7/28. There are additional reasons to choose a fit with the p — > / (600)7 
decay: 

• The branching fraction of the 00 — > 7r°7r°7 decay obtained in the fit III, 
B(lo -> 7T 7r°7) = (11.81^ ± 1.4) x 10~ 5 is above the GAMS result (7.4 ± 
2.5) x 10~ 5 [19] by two standard deviations. The latter result obtained in 
n~p collisions has no p contribution. 

• The value of the coupling constant g puJ7T = 18.6 ± 1.1 obtained in the fit III 
is above that from our analysis above 1 GeV [22] by almost two standard 
deviations. 

• The recent analysis of the process e + e~ — > 7r°7r°7 by SND [21] also showed 
evidence for the p° — > f (600)7 decay. 

Thus, we choose the first model as our final result. The p° — > f (600)7 branch- 
ing fraction is calculated from the results listed in Table 2 taking into account 
that B(/ o (600) -> 7r 7T°) = 1/3. 

3.3 Invariant mass spectrum 

Figure 3 shows the spectrum of a 7r°7r° invariant mass for experimental events 
from the u meson energy range (770-800 MeV). The experimental distribution 
agrees well with the 00 — > p7r° decay model, however, the contribution from the 
uo — > /o (600)7 decay is also acceptable. The existing statistics is not enough to 
distinguish these contributions, therefore we obtain only the total branching 
fraction of the tu — > 7r°7r°7 decay. 

3.4 Search for the decay u — > rjir ^ 

For a search of events of the process e + e~ — > ??7r 7 we first apply the same 
criteria as for the preliminary selection of e + e~ — > 7r°7r°7 events. After that 
a kinematic fit requiring energy-momentum conservation is performed with 
the additional reconstruction of one soft 7r° meson and a good reconstruction 
quality, x 2 < 6, is required. To reject the dominant background from the 
process e + e~ — > 7r°7r°7, we perform an additional kinematic fit with the 7r°7r°7 
hypothesis and reject events that are consistent with it, x^o 7r o 7 < 6. Then 
we look for a possible r\ signal in the invariant mass of two hard photons 
of the remaining three, M 77 . The M 77 distribution is approximated with a 
Gaussian for the signal and a polynomial function for the background. The 
Gaussian mean value and width are fixed from the MC simulation of the signal 
events. The background shape is obtained using the 7r°7r°7 MC. In all energy 
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Fig. 3. The invariant mass of two n mesons in the u> meson energy range. The points 
with error bars represent the experimental data, the solid histogram shows the MC 
simulation of the uj — > p7r° decay, the dashed one corresponds to the to — > /o (600)7 
decay. The hatched histogram is the estimated background contribution. 




Fig. 4. The M 77 distribution for the r]Tr°j candidates. The histogram represents 
experimental events and the solid curve shows the fit result. The dashed curve 
corresponds to the ry7r°7 MC. 

ranges the resulting r/ir ^ signal is consistent with zero. Figure 4 shows the 
M 77 distribution for events from the u resonance region: 381 MeV < -Ebeam < 
401 MeV. The 90% CL upper limit for the number of f]7r 0/ j events is obtained: 
N V7T o^ < 2.4. Using the detection efficiency of 1.3%, we set the following upper 
limit for the e + e~ — * ?]7r 7 cross section: 

a(e + e- -> rjn ^) < 57 pb , (10) 
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Table 3 

Main sources of systematic errors 



Source Contribution, % 

Selection criteria 8 

Background subtraction 6 

Model uncertainty 5 

Luminosity 2 

Trigger efficiency 2 

Radiative corrections 1 

Total 12 

and for the branching fraction of the u meson: 

B(lo -> 777^7) < 3.3 x 1(T 5 . (11) 



3.5 Systematic errors 

The main sources of systematic uncertainties in the cross section determi- 
nation are listed in Table 3. The systematic error due to selection criteria 
is obtained by varying the photon energy threshold, total energy deposition, 
total momentum, and \ 2 ■ The model uncertainty corresponds to different de- 
tection efficiencies for the pn° and /o(600)7 intermediate states. It also 
includes dependence on the /o(600) mass and width. The uncertainty in the 
determination of the integrated luminosity comes from the selection criteria 
of Bhabha events, radiative corrections and calibrations of DC and Csl. The 
error of the NT efficiency was estimated by trying various fitting functions for 
the energy dependence and variations of the cluster threshold. The uncertainty 
of the radiative corrections comes from the dependence on the emitted photon 
energy and the accuracy of the theoretical formulae. The resulting systematic 
uncertainty of the cross section in Table 1 as well as of the branching fractions 
in Table 2 is 12%. 



4 Discussion 

The obtained values of the branching fractions, £>(p° — > n°n ^) = (5.21^3 ± 
0.6) x 1CT 5 and B(lo — > 7r°7r°7) = (6.4l2;o±0.8) x 1CT 5 , are in good agreement 
with the previous measurements by GAMS [19] and SND [21]. Both values 
are significantly higher than those predicted by the vector dominance model 



11 



Table 4 

Predictions for branching fractions of p, lo — ► 7r°7r°7, rjir ^ decays. 

Mode Branching fraction 

(1.1-4.7) xl(T 5 
(1.4-8.2) xl(T 5 
2xl0- 10 -4xl0- 6 
8.3xl0~ 8 -6.3xl0~ 6 

with the standard value of the coupling constant: ~ 1 x 10~ 5 and ~ 3 x 
1CT 5 [4], respectively. An attempt to explain the obtained branching ratios 
results in a high value of g puJ7r contradicting the other observations like, e.g., 
the experimental values of the uj — > 7r + 7r~7r° and uj — > 7r°7 widths. 

Theoretical papers on the subject [1,2,3,4,5,6,7,8,9,10,11,12,13] offer a broad 
choice of effects influencing the discussed decays. In addition to the ujpix 
transition they include: p — uj mixing, pion and kaon loops, various scalar 
(/o(600),/o(980),ao(980)) and tensor (/ 2 (1270) and a 2 (1320)) intermediate 
mesons decaying into nn (rpr). 

Predictions of these models differ rather strongly from each other reflecting 
various approaches applied by their authors [10]. While most of the recent pa- 
pers agree that the observed value of the branching fraction for the p° — > 7r°7r°7 
decay can be ascribed to the intermediate / (600) state, the situation with 
the corresponding lo decay remains controversial. The corresponding ranges of 
the predicted values of branching fractions are summarized in Table 4. Note 
that from the upper limits for the non-cj7r 7r°7r°7 cross section obtained by us 
at higher energy in Ref. [22] a significant contribution from the fo (980)7 or 
/2(1270)7 mechanisms appears not very likely. 

Much higher data samples of the p and lo decays expected in experiments at 
the upgraded collider VEPP-2000 in Novosibirsk [28] will help to significantly 
improve our understanding of their radiative decays. 

From the obtained results on the cross section of the radiative processes 
e + e~ — > 7r°7r°7, r]7T 0/ ~f one can estimate a possible contribution of the pre- 
viously unstudied radiative processes to the leading order hadronic correction 
to the muon anomalous magnetic moment. To this end we first calculate the 
contribution of the process e + e~ — > 7r°7r°7 using the "bare" cross section, a, 
from Table 1 in the energy range below 920 MeV. The result contains a piece 
coming from the uj — > 7r°7r°7 decay already taken into account in Ref. [29] 
in the whole lo meson contribution, = (37.96 ± 1.07) • 10~ 10 . This lo me- 
son contribution is subtracted from the value above using the branching ratio 
B(uj -> 7r°7r°7) with the result (6.08 ± 0.82) • 10~ 12 . A possible contribution 



p o ^ vr o 7r o 7 

00 — > 7T 7r°7 
p° — ► r/7r°7 
lo — ► mr°7 
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from the process e + e — > ty + tt 7 is twice that of e + e — > 7r 7r°7, so that 

aj° l7r7r7 (600 MeV - 920 MeV) = (18.2 ± 2.5) • 1(T 12 . 
Adding the contribution from the r/nj final state, we finally obtain: 

a LO,rad( 60Q Me y _ g 20 Me y) < QM . 1Q -10 ftt g Q% QL 

Adding the upper limit from the energy range 920 MeV - 2000 MeV obtained 
previously [22], we obtain 

a LO,rad( 6Q Q Me y _ 2000 Me y) < o,7 . 1Q -10 ftt g Q% CL 

or about 10% of the current uncertainty of a *f ),had [29]. 



5 Conclusions 



The following results are obtained in this work: 

• Using a data sample corresponding to integrated luminosity of 7.7 pb _1 , the 
cross section of the process e + e~ — > 7r°7r°7 has been measured in the cm. 
energy range 600-970 MeV. The values of the cross section are consistent 
with those obtained by the SND detector [21] and have the similar accuracy. 
The following branching ratios have been determined: B(p — > 7r°7r°7) = 
(5.2^;1 ± 0.6) x 10~ 5 and B{u -> 7r 7r 7 ) = (6.4^ ± 0.8) x IO- 5 . 

• We confirm evidence for the p — > fo (600)7 decay with the branching fraction 
B(P /o(600)7) = (6.0l|f ± 0.9) x 10" 5 reported by the SND Collabora- 
tion [21]. 

• A first search for the process e + e _ — > r]7i°j was performed allowing to set 
the 90% CL upper limits: cr(e + e~ — > i]7i 0r y) < 57 pb in the cm. energy range 
685-920 MeV and B(cu -> ^7) < 3.3 x 10" 5 . 

• A possible contribution of the studied radiative processes to the muon 
anomalous magnetic moment was estimated to be negligible. 

The authors are grateful to the staff of VEPP-2M for the excellent performance 
of the collider, and to all engineers and technicians who participated in the 
design, commissioning and operation of CMD-2. This work is supported in part 
by grants DOE DEFG0291ER40646, NSF PHY-9722600, NSF PHY-0100468, 
INTAS 96-0624, RFBR-98-02-1 117851 and RFBR-03-02-16843. 
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